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We present a complete set of local unitary invariants for generic multi-qubit systems which gives 
necessary and sufficient conditions for two states being local unitary equivalent. These invariants 
are canonical polynomial functions in terms of the generalized Bloch representation of the quantum 
states. In particular, we prove that there are at most 12 polynomial local unitary invariants for 
two-qubit states and at most 90 polynomials for three-qubit states. Comparison with Makhlin’s 18 
local unitary invariants is given for two-quibit systems. 
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Local unitary equivalence is a foundational concept in 
quantum entanglement and quantum information, as it 
provides the key symmetry in classifying quantum entan¬ 
gled states of physical systems [I|. Two quantum states 
are of the same nature in implementing quantum infor¬ 
mation processing if they are equivalent under a local 
unitary (LU) transformation, and maiw crucial proper¬ 
ties such as the degree of entanglement @,01! maximal vi¬ 
olations of Bell inequalities |§-[ll]| , and the teleportation 
fidelity [l3| remain invariant under LU transforma¬ 
tions. Moreover, quantum entanglement in multipartite 
qubits has also figured prominently in many quantum in¬ 
formation processing such as one-way quantum comput¬ 
ing, quantum error correction and quantum secret shar¬ 
ing [2 t 3- this reason, it has been a key problem to 
find a complete and operational procedure to distinguish 
two quantum states under LU transformations. 


In [2, Makhlin presented a complete set of 18 polyno¬ 
mial LU invariants for classifying two-qubit states. There 
are numerous results on LU invariants for three qubits 
states [l^ , some general mixed states pTl - [l^ . , tripar¬ 
tite pure and mixed states [1^. A theoretical method to 
reduce the problem to pure n-qubit states was proposed 
in l2ll| . and later generalized to arbitrary dimensions in 
[2^. From a different viewpoint, 23| gave a procedure 
to find the LU operator for multi-qubits using the core 
tensor method. Very recently a method to judge LU 
equivalence for multi-qubits [26j was also proposed and 
more generally SLOCC inyariants for multi-partite states 
are found [27j| . Nevertheless, it remains a wild problem 
to find a complete set of invariants to answer the LU 
question except for two qubit cases. Even for two partite 
cases it is also desirable to find an alternative set of in¬ 
variants to judge LU equivalence, as the original Makhlin 
invariants contain some nontrivial tensor vectors. 


In this article, we propose a brand new method to 
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quantify polynomial LU invariants for multi-qubit sys¬ 
tems in an operational way. For the special case of two- 
qubit systems, our method is more efficient and needs 
fewer invariants than that in [l4|in general. In fact, we 
show that many invariants in [14| are consequences of 
other invariants, and there are at most 12 invariants to 
determine the LU equivalence for two-qubit states. We 
prove for the first time that there are at most 90 invari¬ 
ants for generic mixed 3-qubit states. We also propose an 
operational method to derive a list of polynomial invari¬ 
ants for generic multi-qubit states. We remark that the 
invariants can not be derived from as the latter aimed 
to compute the LU operator for two equivalent multi¬ 
qubits, while our current work takes a different strategy 
to seek a complete set of polynomial invariants. 

We start our discussion to express an V-qubit state p 
in terms of Pauli matrices a a, a = 1,2,3, 
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where I is the 2x2 identity matrix, crA = J (g) / { 
(Taj, (g) / (g) • • • (g) 7 with (Tq,^ at the jk-lh position and 
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are real coefficients. In particular, Tj = {Tj,Tj, T®), j = 
are three dimensional vectors, Tjk = 

1 < j < k < N, are 3x3 matrices. Generally, Uijs - iM = 
are tensors. 
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Let p and p' be two TV-qubit mixed states. They are 
called local unitary equivalent if 

p' = {Ui ® ... ® Un)p{Ui 0 ... (g) C/at)^ (3) 

for some unitary operators Ui G SU(2), i = 1,2,..., N, 
where f denotes transpose and conjugate. 

Lemma 1 Two mixed states p and p' are local unitary 
equivalent if and only if there are special orthogonal ma¬ 
trices Oi, • • ■ ,On G so(3) such that 

(o,, o,jr,(4) 

for any I < ji < ■ ■ ■ < ju < N, k = 1, 2,..., N. 

Proof. The group SU{2) acts on the real vector space 
spanned by ai, i = 1,2,3 via (^ : 

3 

Okiui, (5) 

where O = {Oki) belongs to 50(3). From ([T]), ([3]) and 

dil) one gets the tensor relation (|1]). Note that this action 
realizes the well-known double-covering map SU(2) -G 
50(3). The sufficiency then follows from the fact that 
50(2) is the universal double covering of 50(3). ■ 

Two-qubit states: To derive explicitly the invariants 
under the transformation ([3|), we first consider the two- 
qubit case. From 0 a two-qubit state is given by the 
3-dimensional real column vectors Ti, T 2 , and the real 
3 X 3-matrix T 12 . Two states p and p' are local unitary 
equivalent if and if there are 50(3) operators Oi and O 2 
such that 

- OiTi, T 2 -O 2 T 2 , 

Ti2 = (0i®02)ri2 = 0iri20*, 

where t denotes the transpose of a matrix. 

We introduce the following sets of 3-dimensional real 
column vectors: 

(Ol) = {Ti, T12T2, ri2T(2Ti, Ti2T(2Ti2T2, • ■ • } C R^ 
(O2) = {T2, Tl^Ti,Tl^T^2T2, Tl^T^2Tl2Ti, • ■ • } C R^ 

which are respectively generated by the (Ti 2 T( 2 )-orbit of 
{Ti,ri 2 T' 2 } and the (T;( 2 '^i 2 )-orbit of {T 2 ,Tl 2 Ti}. Here 
{g) denotes the cyclic group generated by g. By the 
Cayley-Hamilton theorem the minimal polynomials of 
T 12 TI 2 and TI 2 T 12 have degree < 3, therefore it is enough 
to use elements in the orbits up to the quadratic powers. 
It is straightforward to verify that all the vectors in {Of) 
are transformed to Oi{Oi) under the transformation ([b]). 
while all the vectors in {O 2 ) are transformed into 02 ( 02 )- 
Moreover, there are at most three linear independent vec¬ 
tors in {Oi), dim{Oi) < 3, i = 1,2. We say that a two- 
qubit state is generic if dim{Oi) = dim{02) = 3. For 
simplicity, we only deal with generic cases in the follow¬ 
ing. The non-generic (degenerate) cases can be studied 
in details too, see remarks after the proof of Theorem 1. 


Let {pi, ■ ■ ■ , pe} and {vi, ■ ■ ■ , vq\ denote the first six 
(spanning) vectors in {Ox) and (O 2 ), respectively. We 
first give a general result using all the spanning vectors. 

Theorem 1 Two generic two-qubit states are local uni¬ 
tary equivalent if and only if they have the same values 
of the following invariant polynomials: 

{pt,Pj), {vi,vf), z < j = 1,2, • • • ,6. 

. (7) 

tr{Tx 2 Tl 2 r, a = 1,2,3. 

Proof. By using the relations in ([B]), it is direct to verify 
that the quantities given in ([7]) are invariants under local 
unitary transformations. 

For generic states, the matrix T 12 TI 2 is nonsingular, so 
is TI 2 T 12 by trace property. Thus T 12 and are also 

nonsingular. We notice that {Of) {O 2 ) and {O 2 ) 

(Ol) as subsets or subspaces, therefore {Of) ~ {O 2 ) = 
R^ for generic states, and a basis of R^ can be pared 
down from the vectors of (Oi) or {O 2 ) by assumption. 

Assuming that two generic two-qubit states p and p' 
have the same values of the invariant polynomials 0, 
namely, the inner products of any two vectors in {Oi) 
are invariant under p ^ p', one has that there must exist 
an orthogonal matrix Oi such that 

0.{0f) = {0[). 

In particular OiTi = T/. Then we can build the following 
commutative diagram: 

{Ol) ^ {O'x) 


{O 2 ) ^ { 0 ' 2 ) 

Consequently Ti\Oi = O2TI2 in End{'R?), or T{2 = 
0iTi 202- Therefore, p and p' are local unitary equiv¬ 
alent. ■ 

Remark. In the above discussions we are only con¬ 
cerned with the generic case. For degenerate cases, one 
needs to analyze case by case. For instance, let us con¬ 
sider the case Ti = T2 = 0. In this case, dim{Oi) = 0, 
i = 1,2. The only invariants left are tr{Ti2Ti2)°‘, 
0=1,2,3. Note that 

3 

Po=MTi2Ti‘2)“=^A“ (8) 

i=l 

is the oth-power sum of the eigenvalues of Ti 2 T( 2 . A 
well-known result of symmetric polynomials implies that 
any pa {a > 4) is an algebraic function of pi,P 2 , and p^. 
For example, p 4 = ^pf — p\p 2 + \p\ + |piP 3 - Hence 
tr{Ti 2 Tl 2 )°' are invariants for any o > 1. By [ 2 ^ if 
two states p and p' have the same values of tr{Ti 2 Tl 2 )°‘, 
there exists a unitary matrix U such that TI 2 T 12 = 
UTi 2 Tl 2 U\ which means that T 12 TI 2 and TI 2 T 12 have 
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identical eigenvalues. Both T 12 TI 2 and T[ 2 Ti 2 are simi¬ 
lar to diag{Xi, X 2 , Xs). Then there exists an Oi £ SO{3) 
such that T{ 2 ' 7 i 2 = 0iTi2Tl20{. Similarly there ex¬ 
ists O 2 such that = O 2 TI 2 T 12 OI. Subsequently 

^12 = O 1 T 12 OI for some Oi and O 2 , so p and p' are local 
unitary equivalent. 

We now sharpen the result of Theorem [TJ Since there 
are at most three linearly independent 3-dimensional vec¬ 
tors of Pi and Vi in o respectively, one can apply Theo- 
rem[T]to the basis vectors. The standard Gaussian elimi¬ 
nation on the matrix [/zi, • • • , pe] can pare down the col¬ 
umn vectors into a basis {pi-^, pi^, Pi^} of (Oi), where 
{*i!* 2 :* 3 } C {1,2,-■■ , 6 }. This means that the number 
of independent invariants that are used to judge the local 
unitary equivalence of two generic two-qubit states is at 
most 15 in general (instead of 33 as is Theorem [IJ. In 
fact, further analysis can reduce the number to at most 
12 polynomial invariants. 

Theorem 2 Two generic two-qubit states are local uni¬ 
tary equivalent if and only if they have the same values 
for the following 12 invariants: 

{T„iT,2Tl2fT,), {T2,{Tl2T^2fT2), ( 9 ) 

{T^,{T^2Tl2fW2), /3 = 0,1,2, (10) 

tr(Ti2Ti*2)“, a = 1,2,3. (11) 

Proof. The set {Of) is a union of two orbits 
{T 12 TI 2 ) ■ Ti and {T 12 TI 2 ) ■ Tj 2 T 2 - The independent in¬ 
ner products given in Theorem [T] are (Ti, (T‘ 2 Tj 2 )^T'i), 
(Ti, {Tx 2 Tl 2 YTi 2 T 2 ) for ^ = 0,1, 2, 3, due to the Cayley- 
Hamilton theorem and the fact that {T 12 U, v) = {u, TI 2 V) 
for any vectors u,v (T 12 ia a real matrix). Similarly the 
orbit {O 2 ) will only contribute the remaining indepen¬ 
dent inner products (T 2 , {TI 2 T 12 YT 2 ), /3 = 0,1, 2, 3. 

We claim that the 3 invariants with [3 = 3 are not 
needed if the traces m are known. The Cayley- 
Hamilton theorem says that 

(Ti 2T‘2)' = ei(Ti2r‘2)" - e2(ri2Ti*2) + 63/, (12) 

where are the elementary symmetric polynomials in 
the eigenvalues A^. By the fundamental theorem of sym¬ 
metric polynomials, the can be expressed as classical 
polynomials in the traces Pa ( 0 ), i.e. Ci are classical in¬ 
variant polynomials of the density matrix: 

ei=pi, e 2 = i(pi-p 2 ), (13) 

63 = I(Pi- 3P2P1 + 2p3). (14) 

o 


Plugging (fT^ into the three invariants (Ti, (Tf 2 T’i 2 )^Ti) 
etc., we see that they are given by linear combinations of 
the invariants (I MTOll with fixed coefficients of the classical 
invariant polynomials (nm of the density matrix, so 
they are redundant. ■ 

As we commented above if we use the Gaussian elim¬ 
ination we also worry about just 12 invariants, i.e., if 
we add /3 = 3 in the first set of invariants (Emni) for the 
3 basis vectors we can waive the trace identities. Hence 
the total number of invariants is at most 12 either way. 
We still include the trace identities m for the sake of 
general (non-generic) cases. 

Multi-qubit case: To simplify presentation, we intro¬ 
duce the following notation: Ty = Tb. We say that a 
word of Ti,Tij is admissible if the adjacent subindices 
match. For example, Tj 2 r 2 , T 12 T 21 T 1 T 12 are admissible 
ones. 

We first consider the three-qubit case to present our 
general results. In this case, corresponding to o, a quan¬ 
tum state has the form: 


p= -/ ++ VTy +r^23CT(iV(2)a(3). 

o 

1<J 

(15) 

If two states p and p' are local unitary equivalent, then 
there are orthogonal matrices Oi £ SO{3) such that 


T(=Oiri, T^ = 02T2, U = 03 T 3 , (16) 

TI 2 = O 1 T 32 OI T (3 = O 1 T 32 OI T !23 = O 2 T 23 OI (17) 
^123 = (Ol Gl O 2 Gl 03 )Ti23. (18) 

It is known [I^ that the last relation m is equivalent 
to either of the following two relations: 


U 23 — 0iTi23(02G03)*, ^123 ~ ( 01 ( 8 ) 02 ) 712303 . (19) 


Here T 123 is understood as the bipartition Ti |23 (resp. 
7 "i 2 | 3 ) in the first (resp. 2nd) equation of (fTOl) . To state 
our results we introduce two subsets of vectors: 


//o \ _ (rp rp (rp rpt \Brp rp rpt rp rp rpt rp (rp rpt \8 rp / rp rpt \2rp \ ^ Xj 3 

\^lh|23 — 1-^15-^123^23-^ 23J -^23,-^123^123-^15-^123-^ 123-^1231-^23-^ 23/ -123, U 123-1123 J -ll,'”/^^, 
(02 G 03)l|23 = {^23, Tl23TlXl2sTl2s{T2sT^2sfT23,Tl23Tl2sTl23Tl, iT23T*3fT23, • • • } C R® ~ G R^ 
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where P = 0, • • ■ , 3, which are respectively the 

m 23 T* 23 )-Orbit of {Ti,Ti23(r23r2‘3)/5r23|/3 = 0, 1, 2, 3} 
and the {Tl^^Ti2z)-OTUt of {(T 23 rj 3 )^T 23 , r*23Ti| ^ = 
0,1,2,3}. Here r23 is taken as its (column) vector re¬ 
alignment in R® and T '123 is folded as a 3 x 9 -matrix, by 
viewing T123 as the bipartition 1|23 and TI23 is the trans¬ 
pose with respect to such partition. As before we also use 
the same symbols for the corresponding real subspaces. 

Similarly, by permuting the indices we define {O2) '■= 
((^2)2131 and (03) := (03)3112 to be the (1231^231)“ 
orbit of {T 2 ,T 23 i(T 3 ir 3 ‘i)' 5 r 3 i |/3 = 0 , 1 , 2 , 3 } and the 
( 7 ^mT’i 23 )-orbit of {T 3 ,T 3 i 2 (Ti 2 T‘ 2 )'^ri 2 | /3 = 0,1, 2, 3} 
respectively. Here the 3 x 9 -matrix T231 (resp. T312) 
is the realignment of ri23 with respect to the parti¬ 
tion of {123} into {2|31} (resp. {3|12}). Let (Oi) = 
{M1) 1*2,1*3}, (O2) = {* 2 i, 1^2, **3} and (O3) = {Ai, A2, A3}; 
{O2 ® 03)1123 = joi, < 32 , ■ • ■ , Og}, (O3 ® Oi) 2|31 = 

1 ^ 1 ,/ 32 ,..., ^9}, and (Oi (g) 02)3112 = {71,72, •,79}- 


Theorem 3 A three-qubit state p is local unitary equiv¬ 
alent to a three-qubit state p' if and only if the respective 
invariant polynomials are equal: 


{p-ii Pj) — {PiTpj)^ 

(A,A,) = (A',A'), 

(«fc,a/) = {a'k.a'i), 


1 < * < j < 3 
l<k<l<9 
l<k<l<9 
l<k<l<9 


( 20 ) 


Proof. By the result of two-qubit case, the invariance 
of inner products of vectors in (Oi) implies the existence 
of orthogonal matrices Oi, i = 1, 2, 3 such that Eqs. IM 
ini) hold. Thus we are left to show that the orthogonal 
matrices Oi also satisfy Eq. (fT^ or equivalently Eq. (IT^ . 

We use a similar method of Theorem [T] to show this 
by viewing the three-qubit state p as a bi-partite one 
on (g) C® and partition the hyper-matrix T 123 as a 
rectangular matrix T 1123 . Then the 3 x 9-matrix T 123 
maps the subset (O 2 (g 03)1123 into the subset (Oi)i |23 
by left multiplication. 

We have already seen that there exists an orthogonal 
matrix Oi such that 


The following result shows that there are at most 90 
invariants to judge LU equivalence for two three-qubit 
states. 


Theorem 4 Two generic three-qubit states are local uni¬ 
tary equivalent if and only if they have the same values 
of the following invariants: 

(Ti, (Ti 2 T{ 2 )“ri), (T 2 , (ri‘ 2 Ti 2 )“T 2 ), 

(ri,(Ti2Ti‘2)“Ti2T2), 

tr(Ti2Ti*2)^ tr(Ti3r{3)^ triT23T*sf, 

(iri,(Ti|23T3%3)'=Ti), (r2,(r2|3iT2‘33)'=r2), (21) 

{Ts, (ir3|i2T3*|i2)"r3); (T 23 , (r{|23ri|23)'=T23), 

/rri fT't rp \krr\t rr\ \ 

\^23, Ui|23^123j -^1|23^1^ 

^*■(i^l|23’^l|23)^ ^*■(T'2|3l’^2|3l)^ ^*’('73112^^3112)*- 

where a = 0,1,2; /3 = 1,2,3 and k = 0,1, ■■■ ,8; I = 

I,--- ,9. 


The above criteria can be generalized to multi-qubits. 
Define (Oi) = (Ti, T„...„i...i_i|iT„, ■ • •) c as the 
(T,*T.)-orbit, where i = l---i---n means the index 

\ II 11/ ^ 

i is absent. In general for any strict sequence i = 
(H-'Ofc) (i.e. distinct ij’s), we define the (T‘Tj.)-orbit 
(Oi) = (Tj, • • •), where the admissible generating words 
have only i when crossing out redundant strings, e.g., 
T 312 T 12 T 1 is a word of indices 3,1 when crossing out 
12. Then we have the following result. Let (Oi) = 
{Mi, 1 * 2 , -,1*771}, (O2) = {j2i,Z/2,-.-,*^m},---, {On) = 

{Ai, A 2 ,. -., Am}, and more generally, for any strict se¬ 
quence i, let (Oi) = {ri, - - - ,r„}, where n = n(i). Let’s 
list these (Oi) as (Oi^.) = • - , r,^]}, j = 1, • • • , M. 

Theorem 5 Two generic multi-qubit states p and p' are 
local unitary equivalent if and only if the respective in¬ 
variant polynomials are equal: 


{^i 


(1) ..(1) 


\‘i 1 s 


) = {n 


( 1 )' 


r(')') - 




( 22 ) 


where each pair of indices {i,j) are such that 1 < i,j < 
m(i) for the sequences ii, - - - ,im- 


0.(0.) = (O'). 

and Eqs. (El) hold. Then we can directly verify that the 
following diagram is commutative: 

(02 0 03 ) ( 0^0 0 (,) 


(0'i)i|23 

Consequently T 323 (O 20 O 3 ) = O 1 T 123 in En(i(R^ 0 R^), 
o*' '^123 — OiTi23(02 0 03 )*. ■ 



Proof. We use induction on n to reduce the problem to 
(n — l)-partite qubits. Note that for any sequence i of 
indices for n-partite state, we can view the elements in 
(Oi) as (Oi/ 0Oj) where i' is obtained by realignment of 
the Block matrix with respect to the index j, and i' is 
obtained from i after the realignment. Then we can use 
the similar commutative diagram 

(O 20 O 3 ... 77 ) " (O' 0 O'...„) 



(Ol)l|2...77 -^ (0{)l|2...77 
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to get 'ri2...„(02 0 Os - n) = OiTi 2 ...n in EndCR^ (g) 
R3("-2)), orTi' 2 ...„ = OiTi 2 ...„( 02 ® 03 ...n)*. Here O^.-.n 
is an orthogonal matrix in the bigger orthogonal group. 
Then we use the induction to argue further for the ma¬ 
trix Ti| 2 ...„ viewed as a reduced matrix for (n — l)-partite 
state to get the final result. ■ 

Conclusions and Remarks: It is a basic and funda¬ 
mental question to classify quantum states under local 
unitary operations. The problem has been figured out 
in [2lL l22j| for pure multipartite quantum states. How¬ 
ever, it is much more difficult to classify mixed quantum 
states under LU transformations. Operational methods 
have been presented only for non-degenerate bipartite 
states. Although the authors in have shown that the 
problem of mixed states can be reduced to one of pure 
states in terms of the purification of mixed states mathe¬ 
matically, the protocol is far from being operational. We 
have provided an operational way to verify and classify 
quantum states by using the generalized Bloch represen¬ 
tation in terms of the generators of SU{2). We remark 
that [ 2 ^ gives a practical procedure to compute the LU 


operator for two equivalent multi-qubits, but it can not 
derive the polynomial invariants from the procedure, as 
it is based on a different strategy. In our current ap¬ 
proach we set our goal to write down a set of simple 
invariants with which two states can be easily checked 
if they are LU equivalent. Since the coefficients (ten¬ 
sors) in the representation can be determined directly by 
measuring some local quantum mechanical observables- 
Pauli operators, the method is experimentally feasible. 
Our criterion is both sufficient and necessary for generic 
multi-qubit quantum systems, thus gives rise to a com¬ 
plete classification of multi-qubit generic quantum states 
under LU transformations. 
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